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Abstract 

We study projective-anticipating, projective, and projective-lag synchronization of time-delayed chaotic 
systems on random networks. We relax some limitations of previous work, where projective-anticipating 
and projective-lag synchronization can be achieved only on two coupled chaotic systems. In this paper, 
we can realize projective-anticipating and projective-lag synchronization on complex dynamical networks 
composed by a large number of interconnected components. At the same time, although previous work 
studied projective synchronization on complex dynamical networks, the dynamics of the nodes are coupled 
partially linear chaotic systems. In this paper, the dynamics of the nodes of the complex networks are 
time-delayed chaotic systems without the limitation of the partial-linearity. Based on the Lyapunov stability 
theory, we suggest a generic method to achieve the projective-anticipating, projective, and projective-lag 
synchronization of time-delayed chaotic systems on random dynamical networks and find both the 
existence and sufficient stability conditions. The validity of the proposed method is demonstrated and 
verified by examining specific examples using Ikeda and Mackey-Glass systems on Erdos-Renyi networks. 
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In 1999, Mainieri and Rehacek observed projective synchronization in coupled partially 
linear chaotic systems where the drive and response vectors synchronize up to a constant ra- 
tio a (scaling factor). Complete synchronization and anti-phase synchronization are proved 
to be the special cases of projective synchronization in cases of a = 1 and a = —1, re- 
spectively. This proportional feature can be used to extend binary digital to M-nary digital 
for achieving fast communication. With the development of research on complex systems, 
more and more researchers carried out the study about complex dynamical behaviors on 
networks. While most studies focused on complete and phase synchronization in various 
networks, little attention has been paid to projective synchronization. In this paper, based 
on the Lyapunov stability theory, we theoretically analyze both the existence and sufficient 
stability conditions of the projective-anticipating, projective, and projective-lag synchroniza- 
tion of time-delayed chaotic systems on random networks. Numerical experiments for the 
Ikeda system and the Mackey- Glass system show that the control method works. 



I. INTRODUCTION 



Chaos synchronization has attracted considerably increasing attention and become an active 
area of research for both theoretical interests and practical applications, since seminal work of 
Pecora and Carroll Over the last decade, following the complete synchronization [1], several 
new types of synchronization have been found in interacting chaotic systems, such as general- 
ized synchronization Q2[], phase synchronization |3], anti-phase synchronization [4|], projective 



synchronization |5Q, lag synchronization [6] and anticipating synchronization [|7|]. Complete syn- 
chronization is characterized by the convergence of the two chaotic trajectories, y(t) = x(t). 
It appears only when interacting systems are identical. Generalized synchronization means the 
amplitude of the slave's state variable correlated with that of master's by a generic function, 
y(t) = F(x(t)). Phase synchronization is defined as the entrainment of phases of chaotic os- 
cillators, n& x — m&y = constant (n and m are integers), whereas their amplitude remains 
chaotic and uncorrelated. Projective synchronization is the dynamical behavior in which the am- 
plitude of the master's state variable and that of the slave's synchronizes up to a constant scaling 
factor a (a proportional relation). Complete synchronization and anti-phase synchronization are 
the special cases of the projective synchronization in cases of a = 1 and a = — 1, respectively. 
Lag synchronization means a coincidence of shifted-in-time states of two coupled systems, the 
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state variable of the drive system is delayed by positive r in comparison with that of the driven: 
y(t) — x(t — t), (r > 0). Whereas for anticipating synchronization, the driven system anticipates 
the driver: y(t) — x(t + r), (r > 0). Among those synchronization, projective synchronization 
is one of the most interesting problems, because of its proportion between the synchronized dy- 
namical states. In application to secure communications, this feature can be used to extend binary 
digital to M-nary digital communication [|8|] for achieving fast communication. The early study 
of projective synchronization reported that the projective synchronization was usually observable 
only in the coupled partially linear systems [9]. Following work has extended that to a general 



HO. Recently, Hoang et 



class of chaotic systems without the limitation of partial-linearity [fiol 
al. investigated a new synchronization in time-delay chaotic system, which they called projective- 
anticipating synchronization [12]. Projective-anticipating synchronization is a combination of the 
well-known schemes of projective and anticipating synchronization. That is, the driver synchro- 
nizes with the driven under the anticipating synchronization scheme and the amplitude is correlated 
by a scaling factor a: y(t) = ax(t + r), (r > 0). In the case of projective-lag synchronization 



113J], the amplitude of the master's and slave's state variables is correlated by a scale factor a: 
y(t) = ax{t - t), (r > 0). 

However, most existing work about projective-anticipating and projective-lag synchronization 
focused only on two coupled chaotic systems. Systems consisting of many interconnected subsys- 
tems are ubiquitous in nature and social science which can be described by complex networks [14]. 
Complex networks are usually composed by a large number of interconnected components (nodes) 
representing individuals or organizations and edges mimicking the interaction among them. With 
the development of research on complex systems, more and more researchers carried out the study 
about complex dynamical behaviors on networks. While most studies focused on complete Ill5n 
and phase synchronization Jl^ll in various networks, little attention has been paid to projective 
synchronization. Very recently, Hu et al. Ill7ll studied the projective synchronization on drive- 
response dynamical networks by considering coupled Lorenz chaotic systems. But due to finite 
signal transmission times, switching speeds and memory effects with both single and multiple 
delays are ubiquitous in nature, technology and society [18, JjJ]. Time-delayed systems are also 
interesting because the dimension of their chaotic dynamics can be increased by increasing the 
delay time sufficiently Il20ll . From this point of view, these systems are especially appealing for 
secure communication schemes. In addition, time-delayed system can be considered as a spe- 
cial case of spatiotemporal systems [211]. So it is natural to consider an interesting topic whether 
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we can achieve projective-anticipating, projective, and projective-lag synchronization of time- 
delayed chaotic systems on complex dynamical networks? In this paper, inspired by the above 
discussions, we extends the work on projective-anticipating, projective, and projective-lag syn- 
chronization of two coupled chaotic systems to complex dynamical networks. We attempt to 
achieve projective-anticipating, projective and projective-lag synchronization in a general class 
of time-delayed chaotic systems related to optical bistable or hybrid optical bistable device on 
complex dynamical networks. 

The layout of this paper is as follows. In Sec. II, according to the Lyapunov stability the- 



ory 1220, we theoretically analyze both the existence and sufficient stability conditions of the 
projective-anticipating, projective, and projective-lag synchronization of time-delayed chaotic sys- 
tems on complex dynamical networks. In Sec. Ill, the well-known Ikeda and Mackey-Glass sys- 
tems are considered as the dynamics of single node to prove the validity of the proposed theoretical 
approach in Sec. II, respectively. Finally we end this paper by a short conclusion in last section. 

II. THEORETICAL ANALYSIS USING A GENERIC MODEL 

For convenience, we study different types of projective synchronization of infinite-dimensional 



chaotic systems on Erdos-Renyi (ER) random networks [23]. Consider an ER network consisting 
of 1 + iV nodes which constitute the drive-response dynamical networks [|17|l . In response dynam- 
ical networks, we connect each pair of nodes with probability p, and denote the degree of node i 
with h L . Each individual node is a delay-differential system related to the optical bistable or hybrid 



optical device Il24ll 



T 'x(t) = -I3x{t) + fxf(x(t - t R )) (1) 

where x{t) is the dimensionless output of the system at time t, t R is the delay time of the feedback 
loop, t is the response time of the nonlinear medium, fi is proportional to the intensity of the 



incident light and (3 is the parameter. In Eq. (Q~|), f{x) is a nonl 
the system, e. g., f(x) = 7r[l+2_Bcos(x+x )] for Ikeda model [|25| 
Vallee model Il26n . f(x) = sm 2 (x — x ) for the sine-square model [27], and f(x) = ax/(l + x c ) for 



inear function of x, characterizing 

x) = n[A— sin 2 (x— x )] for 



Mackey-Glass model Il28h . In order to observe the different types of projective synchronization, 



the drive-response dynamical networks are described by the following equations: 

T 'x{t) = -f3x{t) + nxf{x{t - n)), (2a) 
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(t - ) N 
r'Ut) = ~PVi{t) + ^f( m Tl> ) + Hf(x(t - r 2 )) + e ]T G ijVj , (z = 1, 2, . . . , N) (2b) 

a i=i 
where x denotes the drive system, ^ denotes the response dynamical network systems, e is the 
coupling strength, a is a desired scaling factor present for the projective-anticipating, projective, 
and projective-lag synchronization, ti is the feedback delay time in the coupled systems and r 2 is 
the coupling delay time between systems x and yi. The graph topology is encoded in the Laplacian 
G, a symmetric matrix with zero row-sum, 



Gij 



-h, if i = j, 

1, if node i connects node j, (3) 
0, otherwise. 



The eigenvalues Aj of G are real and nonpositive I29L 

Throughout this study, for simplification we use the notation x T = x(t—r) . In the following, we 
will investigate different kinds of projective synchronization. One can find that under conditions 

fi = a Hi - fi 2 , (4) 

Eqs. (|2al) and (l2bl) allow for synchronization manifold 

Hi cix T2 _ T1 . (5) 

We define the error systems as the difference between the Eqs. dial) and (l2b~l) : e, = — ai T2 _ T1 , 
one obtains error dynamics 



Ui QLX T , 2 — T1 

1 N 

— [-f3ei + fx 2 f(— 
t a 



— [-(3ei + H2f{^-) +tif(x n ) -afiif(x T2 ) + e^G ij e i ], (6) 

3=1 



according to Eq.©, contained in Eq.© can be expressed as f(x T2 ). Therefore, combining 

Eq.©, Eq.© can be reduced as 

1 N 

ei = — [-fci + e 5^ G v e i\ ■ ^ 
T i=i 
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Consider a Lyapunov function in the form 



N 



He* 

i=i 

n 



(8) 



according to the Lyapunov stability theory [|22Q, if the function © satisfies the following condi- 
tions 

7(e) > ife ^ 0, 



(9) 



and 



7(e) = ife = 0, 
7(e) < ife ^ 0, 



(10) 



7(e) = ife = 0, 

then e will asymptotically converge to zero as time tends to infinity leading to Hindoo \\yi 



0, i.e., projective-anticipating, projective, and projective-lag synchronization of time- 



delayed chaotic systems can be realized on ER networks. 

Surely, the Lyapunov function ® satisfies the condition ®. For the condition (TTOl) . from Eqs. 
© and ® , we get 

TV 

7 = 2J ^iii 
i=i 

a N N N 

Hehee^ 

i=l i=l j=l 

= eG'e T , 



(ID 



where G' = 4-G — 4-J/v, e = (ei, e 2) ....... , e N ). Because G have zero row-sum matrix, accord- 
ing to the stability theory and the Gerschgorin's disk theorem, we can obtain that the sufficient 
stability condition for the synchronized manifold © can be written as 4 > and -4- > 0. Under 
this sufficient stability and existence conditions ©, with the change of magnitude relation between 
feedback delay time n and coupling delay time r 2 , we can observe different types of projective 
synchronization. Let A r = \r 2 — ti\. If r 2 < T\, the synchronization manifold yi = ax T2 - Tl © 
can be written as y^t) = ax(t + A T ), i.e., we can observe projective-anticipating synchronization. 
If r 2 = 7~i, the synchronization manifold y^ = ax T2 _ T1 © can be written as y^t) = ax(t), i.e., we 
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can observe exact projective synchronization. When r 2 > t%, we can derive Ui{t) = ax(t — A T ), 
that is to say, projective-lag synchronization can be observed. 

Remark 1 The value of scaling factor a has no effect on the error dynamics of the system (Eq. 
CO)) because the values of di are independent of the scaling factor a. So we can arbitrarily direct 
the scaling factor a onto any desired value. 



III. NUMERICAL SIMULATIONS 



In this section, we give two illustrative examples. 

Remark 2 Because the whole response dynamical networks realize complete synchronization, 
we can get the projective-anticipating, projective, and projective-lag synchronization between the 
drive system and the response dynamical networks with the same scaling factor a. 

Remark 3 In the first example and the second one, the ER networks consist of 1 + N = 1001 
nodes. In response dynamical networks systems, connectivity probability is p = 0.01. Because 
the size of ER networks is large, it is difficult to get the numerical results between the drive system 
and the whole response dynamical networks. In simulations, we give the numerical results with 
three arbitrarily selected response systems. If they can achieve projective-anticipating, projective, 
and projective-lag synchronization with the drive system x, the other response systems can do it 
with the same scaling factor a. 

Example 1. Consider the nodes of the drive-response dynamical networks are Ikeda systems 
as follows: 

x{t) = —f3x{t) + fiisin(x T1 ), (12a) 



N 

Vi(t) = -/3yi(t) + /i 2 sm(— ) + /j,sin(x T2 ) + e GyUj, (i = 1, 2, . . . , N) (12b) 

i=i 

with (5 > and // 12 < 0. Here, j3 is the relation coefficient for the drive system x and response 
dynamical networks is proportional to the power of the incident light, the delay T\ is 

the time required for light to make a round trip in the cavity. Ikeda model was introduced to 
investigate the dynamics of an optical bistable resonator, playing an important role in electronical 
and physiological study and is well known for delay-induced chaotic behavior [jvL l3oll . 

In the following, we will confirm that the numerical simulations fully support the analytical 
results presented above. The coupling strength in response dynamical networks is set as e = 0. 1. If 
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FIG. 1: Projective-anticipating synchronization for a = 0.5, (3 = 5,t± = 2,T2 = = — 16, fi2 = 

— 0.2, fi = —7.8: (a) The time series of the driver system x(t) (solid line) and the driven systems yi(t) 
(dotted line); (b) synchronization manifold between x(t + A T ) and y(t), A T = \t2 — n|; (c) the time series 
of the error systems a (solid line) and (dotted line). 




FIG. 2: Exact projective synchronization for a = 1.5,(3 = 5,ti = 3, T2 = 3,/ii = —16,^2 = — 0.4, /j, = 
—23.6: (a) The time series of the driver system x(t) (solid line) and the driven systems yi(t) (dotted line), 
(b) the synchronization between x(t) and yi(t). 
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FIG. 3: Projective-lag synchronization for a = —2.0,(3 = 5, t\ = 1,t 2 = 2, \i\ = — 20, // 2 = — 0.5, ii = 
40.5: (a) The time series of the driver system —x(t) (solid line) and the driven systems yi{t) (dotted line); 
(b) the synchronization between x(t — A T ) and yi(t), A T = r 2 — T\. 

t 2 < T\, one can observe the projective-anticipating synchronization for parameters a = 0.5, (3 = 
5, Ti = 2, r 2 = 1, /ii = —16, /i 2 = —0.2, /i = —7.8, with the anticipating time A T = |r 2 — T\\ = 1 
. In Fig. Uta), we can observe that the driven systems anticipates the driver, at the same time, the 
amplitude of x(t + 1) and yi{t) correlates with each other by y^t) = 0.5x(t + 1). Figure [TJb) 
shows the time-shifted plot of x(t + 1) and yi(t). It is clear from the scales of the coordinate axes 
that the slope of the line is 0.5. In Fig. EJc), we give numerical results of the time series of the 
error systems a = i J2iLi IVii^)— < yi(t) > \ and ej, where < • > indicates the average of 
nodes of response dynamical networks. From it we can see that the whole response dynamical 
networks already realize complete synchronization when projective-anticipating synchronization 
among drive-response complex dynamical networks can be observed. With T\ = r 2 , we can 
observe the exact projective synchronization between Eqs. (I12al) and (I12bl) for parameters a = 
1.5,(3 = 5, Ti = 3,r 2 = 3,/ii = — 16, /i 2 = — 0.4, /x = —23.6. From Fig. Eta), one can find 
that the phase angle between the synchronized trajectories is zero. It reduces to the complete 
synchronization if a = 1.0. From Fig. [2[b), we can see the synchronization between x(t) and 
y%(t) for a = 1.5. For coupling delay t 2 being larger than feedback delay n, we can observe the 
projective-lag synchronization for parameters a = —2.0,(3 — 5,Ti = l,r 2 = 2,/ii = — 20, /i 2 = 
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FIG. 4: Projective-anticipating synchronization for a = 0.5, f3 = 6, t\ = 5,T2 = 2,/ii = — 14, /X2 = 
—0.2,// = —6.8: (a) The time series of the driver system x(t) (solid line) and the driven systems yi(t) 
(dotted line); (b) synchronization manifold between x(t + A T ) and y(t), A T = \t% — T\\. 

—0.5, ix = 40.5. In Fig. |3£a), — x(t) instead of x(t) is convenient for our study. The response 
systems lag the state of the drive system with constant lag time A T = r 2 — T\ = 1 and the 
amplitude of the driver's and driven's state variables is correlated by yi{t) = —2x(t — 1). The 
time-shifted plot of x(t — 1) and y(t) is shown in Fig. Ob). It is also clear from the scales of the 
coordinate axes that the slope of the line is —2.0. 

Example 2. Consider the nodes of the drive-response dynamical networks are Mackey-Glass 
systems as follows: 

±(t) = -Px(t) + inx n /{l + x b Ti ). (13a) 

Initially, it was introduced as a model of blood generation for patients with leukemia. x(t) repre- 
sents the density of circulating cells at time t, when it is produced, and x n is the density when the 
"request" for more blood cells is made. Later, this system became popular in chaos theory as a 
model for producing high-dimensional chaos to test various methods of chaotic time-series analy- 
sis, controlling chaos, etc. The electronic analog of this system has been proposed by Pyragas and 
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FIG. 5: Anti-phase synchronization for a = — 1.0, (3 = 6,ti = 5, T2 = 5, fii = —16,^2 = — 0.2,/i = 
16.2: (a) The time series of the driver system x(t) (solid line) and the driven systems yiit) (dotted line), (b) 
the synchronization between x(t) and yi(t). 




FIG. 6: Projective-lag synchronization for a = 1.5, (3 = 6,ti = 3, T2 = 5,/xi = —16, ^2 = —0.2,// = 
—23.8: (a) The time series of the driver system x(t) (solid line) and the driven systems yi(t) (dotted line); 
(b) the synchronization between x(t — A r ) and A T = T2 — ri. 
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his collaborators [31]. The response dynamical network systems are written as 



N 

Vi(t) = -PUi(t) +/i 2 (?/i, T1 /a)/(l + (y iyT Ja) b ) + /j,x T J(l + x b T2 ) + e^G^-, 

(z = l,2,...,iV). (13b) 



In our simulation, the coupling strength in response dynamical networks is set as e = 0.2. With 
r 2 being less than n, one can observe the projective-anticipating synchronization for parameters 
a = 0.5,(3 — 6, Ti = 5,r 2 = 2,/ii = — 14,/i 2 = —0.2,// = —6.8, with the anticipating time 
A T = |t 2 — T\ I = 3 . In Fig. ID^a), We observe that the driven systems anticipates the driver, at the 
same time, the amplitude of x(t+3) and correlates with each other by yi(t) = 0.5x(t+3). Figure 
@Jb) shows the time-shifted plot of x(t + 3) and y { . It is clear from the scales of the coordinate axes 
that the slope of the line is 0.5. With t 2 = t\, one can observe the exact projective synchronization 
for parameters a = — 1.0, (3 — 6, n = 5, r 2 = 5,//i = — 16, /i 2 = — 0.2, /x = 16.2. When 
a = —1.0, we observe the anti-phase synchronization shown in Fig. Oa), synchronized chaotic 
attractors with an anti-phase pattern where the phase angle between the synchronized trajectories 
is 7r. From Fig. Ob), we see the anti-phase synchronization between x(t) and yi(t). With r 2 being 
larger than r±, one can observe the projective-lag synchronization for parameters a = 1.5,(3 = 
6, 7~i = 3, r 2 = 5, /ii = —16, /i 2 = —0.2, /j = —23.8. In Fig. [(J a X We can observe that in this case 
the response systems lag the state of the drive system with constant lag time A T = r 2 — t\ =2 
and the amplitude of the driver's and driven's state variables is correlated by y^t) = 1.5x(t — 2). 
Figure^b) shows the time-shifted plot of x(t — 2) and yi(t). 

IV. CONCLUSION 

To summarize, we have analytically estimated and numerically simulated projective- 
anticipating, projective, and projective-lag synchronization properties of time-delayed chaotic 
systems related to optical bistable or hybrid optical bistable device on random dynamical net- 
works. Our work may lead to several advantages over existing work:(l) It is capable of realiz- 
ing projective-anticipating and projective-lag synchronization of a general class of time-delayed 
chaotic system related to optical bistable or hybrid optical bistable device on random networks, 



while previous work 11121 . 1 X 311 only mentioned projective-anticipating and projective-lag synchro- 



nization on only two coupled chaotic systems. (2) Although previous work Ill7h has achieved 
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projective synchronization on complex dynamical networks, the dynamics of the nodes are cou- 
pled partially linear chaotic systems; we realize projective synchronization of time-delayed chaotic 
systems on random dynamical networks, without the limitation of the partial-linearity; so it can 
be considered as an extension of the dynamics of each individual node from partially linear 
chaotic systems to non-partially-linear chaotic systems, or an extension from finite-dimension 
to infinite-dimensional chaotic systems. According to the Lyapunov stability theory, we have 
achieved the sufficient stability and necessary conditions for the projective-anticipating, projec- 
tive, and projective-lag synchronization manifolds . We derive that the transition among projective- 
anticipating, projective, and projective-lag synchronizations can be achieved by adjusting the mag- 
nitude relation between the feedback delay time and the coupling delay time. The validity and 
feasibility of our method have been verified by computer simulations of Ikeda and Mackey-Glass 
systems. 
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